Dark solitons in ferromagnetic chains with first- and second-neighbor interactions 
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We study the ferromagnetic spin chain with both first- and second-neighbor interactions. We 
obtained the condition for the appearance and stability of bright and dark solitons for arbitrary 
wave number inside the Brillouin zone. The influence of the second-neighbor interaction and the 
anisotropy on the soliton properties is considered. The scattering of dark solitons from point defects 
in the discrete spin chain is investigated numerically. 
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O ■ I. INTRODUCTION 

<N 

Considerable attention has been devoted to the study of solitary waves in solids. A large amount of theoretical 
and experimental research has been dedicated to one-dimensional magnetic systems [1]. Models corresponding to real 
quasi-one-dimensional magnets are broadly investigated as first they were considered as somewhat simpler than the 
, really interesting three-dimensional systems but then turned out to be interesting in their own right. There was the 
£\j ■ availability of rear magnetic compounds which because of extremely small coupling between neighboring chains, could 
be considered as reasonable realizations of magnetically one-dimensional materials. Soliton solutions for classical 
sine-Gordon chains [2] and Heisenberg chains with various anisotropics [3-6] were obtained and analyzed. In recent 
years there is a renewed interest to the topic due to their application. Spin-wave dark solitons were predicted and 
experimentally generated [7,8]. Solitons in magnetic thin films [9] and in ferromagnets with biquadratic exchange [10] 
were investigated. 

Bright and dark solitons as exact solutions of the nonlinear Schrodinger (NLS) equation have been studied intensively 
[3,4,11]. But, the models describing the physical effects in solids are mainly discrete and the question arises how the 
discreteness modify the properties of the nonlinear excitations [12,13]. Another interesting topic of research with 
practical importance is the interaction of solitons with impurities [14-18]. The influence of the interaction with 
second-neighbors [19] and third-neighbors [20] on the bright soliton dynamics has been also considered. 

In the present paper we study the condition for the appearance of dark solitons in an anisotropic ferromagnetic 
chain with first- and second-neighbor interactions. The propagation for arbitrary wave number is considered and the 
scattering on point defects is investigated. 



II. HAMILTONIAN OF THE SYSTEM 

We consider a ferromagnetic Heisenberg chain of N spins with magnitude S described by the following Hamiltonian: 
- 1 — i i j 



~ T ^(^w $n+2 + Sn+2) ~ ^2 ^ ^n^n+2 ~ ^^(^n) 2 i (1) 
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where both nearest-neighbor and next-nearest-neighbor exchange interactions are taken into account. The spin oper- 
ators S^;, obey the commutation relations 

[S±,S*]=T-S±<%, [Sf,Sr] = 2SfS ij . (2) 

Ji > and Ji > (i = 1, 2) are the exchange integrals and A is the on-site anisotropy constant which can be positive 
(easy axis) or negative (easy plane). If Ji ^ Ji the inter-site anisotropy is included. For Ji = Ji and A = the model 
is isotropic. 

The equations of motion for the operators S£ yield (ft = 1): 



dS+ 
dt 

+ [Ji(Sn+i +Sn-i) + "M^W-a + ^-2)]^ +2AS*S n . (3) 



- -AS+ - 5^[Ji(5+ +1 + 5+_j) + J2(S+ +2 + S+_ 2 )] 
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In the quasiclassical approximation, where th e components of the spin operators S n are complex amplitudes, 
a n = S+/S, a* n = S~ /S and S*/S = \Jl — \a n \ 2 , we have 

i-gjr = -Aa n - [JiS(a n+ i + a„_i) + J2S(a n+2 + a„_ 2 )] \/l - \a n \ 2 

+ a n [JxS (Vl - K+i| 2 + y/l- |^-i| 2 ) + J2S (Vl - K+2I 2 + y/l- |«„- 2 | 2 )] 

+ 2ASa„v / l _ KF- (4) 

The set of differential equations (4) describes our system. In the continuum limit (wide excitations), it transforms 
into a perturbed NLS equation, where the perturbing terms are not only due to the discreteness but also to the 
complicated nonlinear interactions. For J 2 = J\ = J 2 = A = and wide excitations, (4) turns to the discrete 
Ablowitz-Ladik equation which is also completely integrable and has soliton solutions. So, for narrow excitations it is 
necessary to investigate the whole set (4) . Equations similar to (4) has been derived for exciton solitons in molecular 
crystals [21]. 



III. SOLITON SOLUTIONS 



We shall look for solutions in the form of amplitude-modulated waves 

a n (t)=<p n (t)e i ( kn -" t \ (5) 

where k and w are the wave number and the frequency of the carrier wave (the lattice constant equals unity) and 
the envelope <p n (t) is a slowly varying function of the position and time. In the continuum limit, when the soliton 
width is much larger than the lattice spacing (L ^> 1), equation (4) transforms into the following NLS equation for 
the envelope: 



i^ +v ° d £)=( £ -^- bk & + ^, (6) 

where 

t = tS , fl — uj/S, gk — Ji cos A: — Ji + J 2 cos2fc — J 2 — A , s = —A/S — 2gk, 
v g = 2(Ji sin/c + 2J 2 sin2fc) , bk = J\ cos k + 4J 2 cos2fc . (7) 

v g is the group velocity and bk describes the group-velocity dispersion of the carrier waves. Depending on the sign 
of the expression bkgk equation (6) possesses soliton solutions of different types. For negative values bright solitons 
exist while for positive values dark solitons are possible. 

Let us consider the dependence of the coefficients bk and gk on k (Fig. 1 and 2). The group- velocity dispersion bk 
depends only on Ji and J 2 and is not influenced by the anisotropy of the system. Without second neighbor interaction 
(J 2 = 0) bk changes sign at the middle (k = ir/2) of the Brillouin zone. The inclusion of weak second-neighbor 
interaction (J 2 < 1/4) has influence only in the moving of this point to smaller k- values. Strong second-neighbor 
interactions (J 2 > 1/4) lead to a second point where the function bk changes again sign (Fig. 1). The nonlinear 
coefficient g k depends on the anisotropy of the system. For the isotropic case without second neighbor interaction 
(Ji = Ji, J 2 = J 2 = A = 0) gk is negative in the whole Brillouin zone. So in this case equation (6) has bright-soliton 
solutions for < k < n/2 and dark-soliton solutions exist for ir/2 < k < ir. The inclusion of positive on-site anisotropy 
(A > 0) preserves this behavior. The inclusion of negative on-site anisotropy (^4 < 0) leads to the appearance of a wave 
number k c so that gk is positive for < k < k c and negative for k c < k < n [(Fig. 2(a)]. Then the Brillouin zone is 
divided in three regions where the solution alters from dark- to bright- and again to dark-soliton. If k c = ir/2 then the 
solution of equation (6) is of the dark-soliton type in the whole Brillouin zone. The inclusion of weak second-neighbor 
interactions has influence only on the value of k c and the size of the regions where the different solutions exist. Strong 
second- neighbor interactions with J 2 > 1/4 will lead to the appearance of the additional region near the band edge tt 
where the function bk becomes positive (Fig. 1) and consequently the solution is of the bright-soliton type. Fig. 2(b) 
shows the function gk for the case of inter-site anisotropy (Ji = J 2 = A = 0). When the second-neighbor interactions 
are not taken into account the solution of equation (6) is of the dark-soliton type in the whole Brillouin zone. The 
inclusion of second-neighbor interactions will lead again to three or four (if J 2 > 1 /4) regions where different soliton 
solutions exist. 
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FIG. 1: Dependence of the group-velocity dispersion bk on the second-neighbor interaction. Ji = 1; J2 = (solid line), 
J2 = 0.25 (dotted line), J2 = 0.5 (dashed line), and J2 = 1 (dashed-dotted line). 
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ivj. z.. i\uminear coefficient gu for (a) on-site anisotropy (Ji = Ji, A — —0.5) and (b) inter-site anisotropy (Ji — A — 0); 
= 1,2. J2 = (solid lines), J2 = 0.25 (dotted lines), J2 = 0.5 (dashed lines), and J2 = 1 (dashed-dotted lines). 
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In what follows, we shall consider dark-soliton solutions which appear for 

bk9k > 



(8) 



and have nonvanishing boundary conditions, |(^(a:)| 2 — > const at x — > ±00. Note that the condition (8) depends not 
only on the anisotropy constants but also on the wave number k. The dark soliton has the form 



x — vi~ 

ip(x, t) = t^otanh — h iB , 



(9) 



where L and v are its width and velocity. B determines the value of the amplitude at the center, i.e the minimum 
intensity of the soliton |</?(0)| 2 = B 2 , while the maximum intensity is |^(±oo)| 2 = ip 2 + B 2 . As independent parameters 
of the soliton we can choose L, k and B. Then <po, il and v are given by the relations: 
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1b k —— — 

Pv = ~^j2> ^ = £ + 9k(<Po + B 2 ) , v = v g +v d , v d = B^/2g k b k . (10) 

The dark soliton can be considered as an excitation of the background (carrier) wave. Its velocity consists of the 
group velocity v g and a contribution Vd which depends on B. For fixed (<£>q + B 2 ) and k solitons moving faster have 
smaller amplitude ip n . 

The bright soliton solution which exists for b k gk < and |</?(a;)| 2 — > at x — >• ±oo has the form 

<p(x,r) = y sech X VT (11) 



L 



with 



^ = "^' n = £ -J7^ v = v °- (12) 

We like to make some remarks about the influence of the next-nearest-neighbor interaction on the form and velocity 
of the soliton solution. The dependence of the velocity on k is shown on Fig. 3. The group velocity v g is defined 
in the whole Brillouin zone and coincides with the velocity of the bright solitons and of the dark solitons for B = 
[Fig. 3(a)]. It is independent of the anisotropy of the model. The small second neighbor interactions change only the 
value of the velocity, while strong second neighbor interactions (J2 > 1/4) can change also its sign. The velocity of 
the dark soliton itself v d is shown on Fig. 3(b). It is nonequal zero in the regions where the dark solitons exist and 
its sign is determined by B as given in (10). 

We performed numerical simulations based on the general discrete equation (4) and we investigate mainly the dark 
soliton solution. As initial function we put for t = the form 



M0 ) = ( V!|M tanh ^ +lB j e ikn (13) 
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which is exact solution of the quasicontinuum equation, and observe the evolution for different anisotropic cases. For 
wide solitons (L 3> 1) this solution remain stable and its form and velocity are preserved. When the soliton width 
decreases the discreteness effects of equation (4) become important and the parameters of the dark soliton change. 

In what follows we have L — 10 and J\ = 1. Fig. 4 shows the propagation of the dark solitons in the case Ji = A = 
for two fc-values (fci = 0.0157 and fc 2 = tt — ki). When the second neighbor interaction is absent ( J 2 = 0) the solitons 
propagate with the same velocity and amplitude [Fig. 4(a), (a')]. When the second neighbor interactions are present 
the solitons have different amplitudes and velocities [Fig. 4(b), (b') for J 2 = 0.1]. As can be seen from the fc-dependence 
of the velocity from Fig. 3 the value of the velocity for the greater k is smaller. The increase of the second neighbor 
interaction leads to the turning of the solution for fc 2 to a bright soliton type [Fig. 4(c), (c') for J 2 = 0.3]. 

Similar results have been obtained for the case with on-site anisotropy J\ = J%, J% = J2, A ^ 0. 

We like to point out that Eq. (4) is invariant with respect to the transformation a n —> a*e ?7TO when J 2 — > — J 2 , 
Ji — > —J%, J 2 — > —J2, and A — > —A. So the results on Fig. 4(b), (c) can be obtained also for fc 2 provided J 2 has the 
values -0.1 and -0.3, respectively. 



IV. SCATTERING ON POINT DEFECTS 



We studied the interaction of the dark solitons with impurities in the chain whose size is small compared to the 
soliton extent. Then the equation of motion becomes 

= -Aa n - [Ji5(a n+ i + a„_i) + J 2 S*(a„+ 2 + a„_ 2 )]\/l - |a„| 2 

+ a n [,hS (y/i - K +1 | 2 + y/1 - K-il 2 ) + J 2 S (V 1 - K+2| 2 + V 7 ! - Wn-2?)] 

+ 2ASa„ v / l - |a„| 2 + dS n<no a n . (14) 
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FIG. 5: Scattering of a dark soliton for inter-site anisotropy ( Ji = J2 = A — 0), J2 = 0.1 and k\ = 0.0157 [(a), (b)], ki — n — k\ 
[(a'), (b')]- Transmission for d = 0.0152 (a) and d = -0.00652 (a'). Reflection for d = 0.0153 (b) and d = -0.00653 (b'). The 
time is in units of 10 3 / J\S. 

The parameter d describes the strength of the impurity at site no, which can be attractive or repulsive. In our 
numerical simulations we investigate the interaction of a propagating dark soliton of the form (13) with a defect 
placed at the middle of the chain no = 500. We have observed that in the case of attraction the solitons pass through 
the defect for arbitrary initial velocity. More interesting is the interaction with repulsive impurities. 

In the case of repulsion the soliton can be transmitted or reflected by the impurity. For a given initial velocity vq 
there is a critical value of the defect |cf c | below which the soliton is transmitted and above which it is reflected. We 
consider dark solitons which velocities are determined only by k (B = 0). Then, due to the symmetry properties 
of equations (4) and (14) mentioned above, d c > for small k-values and d c < for k near it. Fig. 5 shows the 
scattering patterns for the case of inter-site anisotropy {J\ = J2 = A = 0) and different wave numbers (ki = 0.0157 
and &2 = 7r — fci). When the interaction with the second neighbors is absent ( J 2 — 0) the scattering mechanisms for 
k\ and &2 are equal, as the values for the initial velocity Vo = 0.031 the values for |d c | are the same and it holds: 
0.0108 < \d c \ < 0.0109. When the interaction with the second neighbors is included (J2 = 0.1) the initial velocities 
for the two cases become different and the the values for d c are different, respectively. As the initial velocity for k\ 
(vq = 0.044) is greater than the initial velocity for k% (vq = 0.019) the value for |ef c | is higher. We observe that the 
scattering mechanisms depends on the initial velocity and the strength of the impurity. The anisotropy of the chain 
influences the amplitude of the soliton in this case. 

V. CONCLUSION 

We have investigated the existence and stability of dark solitons in a discrete ferromagnetic chain with both inter-site 
and on-site anisotropy. We have considered the influence of the second-neighbor interactions on the soliton dynamics 
in the whole Brillouin zone. The complicated dependence of the dispersion and the nonlinear coefficients lead to 
regions in the Brillouin zone where strong second-neighbor interactions can turn the type of the soliton solution (from 
bright to dark or vice versa). The scattering of the dark solitons on point defects is also investigated. 
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